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The CW approximation takes into account electrostatic self-interaction contained in the Hartree 
potential through the exchange potential. However, it has been known for a long time that the 
approximation contains self-screening error as evident in the case of the hydrogen atom. When 
applied to the hydrogen atom, the GW approximation does not yield the exact result for the electron 
removal spectra because of the presence of self-screening: the hole left behind is erroneously screened 
by the only electron in the system which is no longer present. We present a scheme to take into 
account self-screening and show that the removal of self-screening is equivalent to including exchange 
diagrams, as far as self-screening is concerned. The scheme is tested on a model hydrogen dimer 
and it is shown that the scheme yields the exact result to second order in {Uq — U{)/2t where f/o 
and Ui are respectively the onsite and offsite Hubbard interaction parameters and t the hopping 
parameter. 



PACS numbers: 71.10.-w, 71.27.+a, 71.15.-m 

I. INTRODUCTION 

In the Hartree approximation^, a system of electrons 
move in a common potential arising from the electro- 
static field of the electrons, in addition to the external 
field. In this approximation, a given electron experi- 
ences the electrostatic potential from the other electrons 
as well as from itself because the common potential or 
the Hartree potential contains the field from the elec- 
tron itself. This unphysical self-interaction is removed 
when exchange interaction is included, leading to the 
Hartree- Fock approximation (HFA)^. In density func- 
tional theory2,^ Perdew and Zunger introduced the con- 
cept of self-interaction correction^ to remove a similar 
problem in the local density approximation (LDA). 

For many-electron systems, such as solids, it is well 
known that the HFA is not satisfactory because it com- 
pletely neglects screening which is very crucial in describ- 
ing the electronic structure of many-electron systems. 
Thus, for example, the Hartree-Fock band gaps of semi- 
conductors and insulators are much too wide and when 
the HFA is applied to metals the density of states at the 
Fermi level becomes unphysically zero due to the loga- 
rithmic singularity in the derivative of the one-particle 
energy with respect to the A;-vector at fc = kp^. The 
simplest known and successful method beyond the HFA 
that cures the band-gap problem and the anomaly of the 
HFA in metals is the GW approximation (GWA)''.!l. The 
GWA includes the effects of frequency-dependent screen- 
ing from first principles and the self-energy in space-time 
representation is approximated by a product of the Green 
function G and the screened interaction W. 

The GWA includes the exchange potential so that it 
is self-interaction free. However, it is contaminated by 
" self-screening" , namely, an electron screens itself, anal- 
ogous to "self-interaction" where an electron interacts 



with its own electrostatic field. This undesirable self- 
screening effect has been a long-standing problem and 
thought to be a source of significant errors in the elec- 
tronic structure. The self-screening problem may be il- 
lustrated by the famous case of the hydrogen atom. Since 
there is only one electron, it is clear that the one-particle 
removal energy or the hole energy is simply given by 13.6 
eV, the Is orbital energy. The Hartree approximation 
applied to the hydrogen atom would yield a too low re- 
moval energy due to the self-interaction error while the 
HFA would give the correct result. Embarassingly, when 
the GWA is applied to the hydrogen atom, it yields a 
wrong result because, as a consequence of self-screening, 
the correlation part of the hole self-energy in the GWA 
is not zerc^. Evidently, since there is only one electron, 
upon removal of the electron there are no other electrons 
that can screen the remaining hole so that the hole self- 
energy ought to be zero. 

The self-screening error is believed to be responsible 
for a number of well-known problems. It has been sus- 
pected for a long time that the presence of self-screening 
in the GW self-energy may be responsible for errors in 
the quasiparticle energies of localized states. It has been 
found that GW quasiparticle energies of core or semicore 
states usually lie above the experimental values. It is ar- 
gued that in the HFA the quasiparticle energies are too 
low due to the absence of screening and when screening 
is taken into account within the GWA, these energies are 
pushed up too high, an indication of overscreening due 
to self-screening. In molecules, a recent comprehensive 
and systematic study of 34 molecules has found that the 
GWA overscreens the Hartree-Fock ionization potential 
leading to underestimation by 0.4 ~ 0.5 eV compared 
to experiment^. In many materials, the energy position 
of the core or semicore states is usually too high in the 
LDA due to self-interaction. GW calculations on the 3d 
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semicore states of a number of semiconductors such as 
GaAs and ZnSe improve the LDA results but the remain- 
ing error is still significantiiS. It is very likely that this 
error owes its origin from self-screening. From physical 
consideration the self-screening error is expected to be 
significant when the states are rather localized but less 
important in extended states 

In this paper, we develop a new scheme which aims at 
correcting the self- screening error in the GW self-energy 
as well as the linear density-density response function 
within the random-phase approximation (RPA)-'^^. An 
interesting consequence of the proposed scheme is the 
fact that the screened interaction W becomes explicitly 
spin dependent, in contrast to the original GWA where 
the screened interaction is spin independent. We also 
furnish a theoretical support for the scheme by showing 
from diagrammatic consideration that the removal of the 
self-screening terms is partially equivalent to adding ex- 
change diagrams. In other words, the self-screening terms 
are cancelled by corresponding terms in the exchange di- 
agrams. 

As an illustration of our scheme, we calculate the 
bonding-antibonding gap of a model hydrogen dimer. We 
have chosen this model because the exact result is known 
allowing for rigorous comparison. Moreover, the calcula- 
tions can be performed analytically so that possible nu- 
merical errors are eliminated and the simplicity of the 
system permits us to analyze the results without unnec- 
essary complicating factors. It is found that the self- 
screening corrected GWA reproduces the exact result to 
order [{Uo — Ui)/2t]^, where Uq and Ui are respectively 
the onsite and offsite Coulomb energies and t is the hop- 
ping integral. 



II. GW APPROXIMATION WITH 
SELF-SCREENING CORRECTION 

A. Theory 



E^(rt,r'i') = iGl{rt,r't')W{r't',rt) 

= iY,gma{rt,r't')W{r't',rt), (4) 

m 

where W is the screened interaction 

W = e-\ (5) 

with e being the dielectric matrix. 

Consider an electron occupying an orbital ipmir prop- 
agating from {r't') to (rt) represented by gma- Another 
electron with the same spin cannot occupy the orbital 
ipma and therefore gma should not participate in the 
screening process during the propagation of the elec- 
tron. Therefore the screened interaction W should be 
calculated using a polarization propagator that does not 
include gma-- However, an electron in the same orbital 
but with opposite spin can naturally participate in the 
screening process. The self-energy then takes the follow- 
ing form: 

S,(ri,r't') = i^5„,(ri,r't')T^„,(rY,ri), (6) 

m 

where 

W.ma =V + vRmaV + , (7) 



The polarization Pma is defined as the polarization with- 
out gma, i-e., no Green function line in P^a contains gma- 
In other words. 



The first step of the procedure is to decompose the non- 
interacting Green function into its orbital components: 



G°(r,r';L.)=^.g„,(r,r': 



gna{r,r';u}) = 



y^na{r)ipt,Ar') 



UJ - £r. 



(1) 



(2) 



Pma — iiGmaGma G*— crG*— (j), (9) 

where Gma is the Green function without gma, namely, 

Gma = Gcr — gma- (10) 

In Fig. [T] the self-energy diagrams corresponding to ([5]) 
are compared with the conventional GW diagrams. 

The correlation part of the GW self-energy with self- 
screening correction is given by 



hi-Pr. 



(3) 



where /i is a one-particle Hamiltonian, Sna ^ £na + 
for an occupied state and £na — >■ £na — for an unoc- 
cupied state. We refer to {gna} as orbital Green func- 
tions. In the GWA the self-energy is given by, using a 
non-interacting G°, 



-^gma{Y, r'; u + cj') W,;(r', r; u') 



27r W + W' - £ma + j<5sgn(£m^ - [l) ' 



(11) 
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FIG. 1: Comparison between the self-energy diagrams in the 
conventional GWA (a) and the GWA with self-screening cor- 
rection (b). In the latter, the screened interaction depends 
both on the orbital and spin of the electron represented by 
(?„, as discussed in the text. 



Writing the correlation part of the screened interaction, 
W, in its spectral representation 



Tire / / /^ a" J ri Dmcr{r' ,r]UJ 
W^A^ ,r;uj)= duj 



duj' 



uj' — uj" — iS 
oj'-oj" + iS ' 



(12) 



the frequency integral over uj' can be performed analyt- 
ically. The correlation part of the self-energy may be 
divided into two parts, 1]°'=" and 



Sr(r,r';c.) 



E 



^^,„ (/?,na(r)An^(r', r; uj")ip;^„{r') ^^^^ 



It is worth noting that the self-screening correction in- 
troduces spin dependence in the screened interaction W 
as can be seen in (jH). Each electron experiences a differ- 
ent screened interaction Wma which is not only orbital 
dependent but also spin dependent according to the or- 
bital occupied by the electron as well as the spin of the 
electron. 

Since in the exact set of Hedin's equations the screened 
interaction W is spin independent the appearance of a 
spin-dependent screened interaction seems unnecessary. 
It is interesting to make comparison with density func- 
tional theory. In principle, the total energy is obtainable 
from the ground-state electron density, which is the sum 
of spin up and down components. In practice, for spin- 
polarized systems it is more favorable to introduce the 
spin variable and regard the total energy as a functional 
of the up and down spin densities. The separation of the 
density into the up and down components mimics the 
true system and captures the essential physics so that a 
relatively simple approximation, such as the local spin 
density approximation, still works well. A presumably 
much more complicated functional would be required to 
achieve the same level of accuracy for the total energy 
if the total density were to be used instead. A sim- 
ilar situation arise in our case, where the orbital and 
spin-dependent screened interactions closely mirror the 
physical situation and thereby promotes a better self- 
energy within the simple GWA. If we kept the conven- 
tional screened interaction, we would need to include ex- 
change diagrams as vertex corrections to cancel the self- 
screening terms, as shown in a later section. It is much 
simpler to remove the self-screening terms than to include 
vertex corrections. 



UJ + uj" 



iS 



I]~(r, r';uj) 



E 



duj 



UJ ~ uj" — e,na + iS 



(14) 



Dma is the spectral function of W^'^jg. and we have used 
the relations 



D^4-iu) = -Dma{^), W^,,,{-UJ) = W^i^A^) (15) 



and 



1 



ImW;';,^(w)sgn(w). 



(16) 



B. Self-screening correction in extended states 



For extended states, the self-screening correction tends 
to vanish. However, from the physical point of view, 
we expect that the self-screening correction is significant 
when the state originates from a localized orbital such 
as the case with the states originating from the 3d or 4f 
orbitals. Consider expanding a given Bloch state in its 
Wannier representatiorti^ 



k„a(r) = exp(ik • R)xr„<t(i-)- (17) 

V ^ R 



The corresponding expressions for the self-energy in the 
conventional GWA are the same as above except that 
Dma is replaced by the spectral function of W instead. 



The Green function is, with ej^na — ^ ^^Jina+iv for occupied 
states and ekno- — >■ Ekno- — iij for unoccupied states, 
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kn 
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TV 



EE 

kn RR' 



exp[ik • (R - R')]xR„^(r)xR.„^(r') 



^ kn R^^R' 



AT EE 



XRna(r)XR„<T(l"') 



kn R 



(18) 



We apply the self-screening correction to the component 
of G° corresponding to R = R'. As before we define 



^ncr — "-^o- yn 



(19) 



5„,(r,r';w) = ^EE 



k R 



XRna(r)x^tna(''0 

^ ~ £kna 



R 



(20) 



where 



(a) exchange 



(b) direct 



FIG. 2: The second-order exhange and direct self-energy and 
polarization diagrams. As shown in the text, the exchange di- 
agrams cancel the self-screening terms in the direct diagrams. 



the corresponding terms in the exchange diagrams in a 
similar fashion as for the first-order self-energy or the 
HFA. We will illustrate the idea for the second-order self- 
energy but it is clear that the argument applies to any 
order. The second-order exchange and direct diagrams 
are shown in the upper part of Fig. [2] 

According to the Feynman rules^"^, the second-order 
exchange self-energy for a given spin is 



S^(a;i,X2) (i)'^ J dx3dx4G{xi,X3)G{x3,X4) 

X G{xi,X2)v{xi - X4)v{x3 - X2), (24) 
where x = (r, t) and 



E 



^ — £kncr 



(21) 



In practice, there may be a problem due to the non- 
analytic behaviour of eknir as a function of k. For very 
narrow band such as the one formed by semicore states 
we may make the following approximation 



N ^ 

R 



XR «<T(r)xRn^(rO 

W - (encr) 



(22) 



v{x-x')^v{v-v')&{t-t'). (25) 
Fourier transformation with respect to t ~ t\ ~ t2 yields 

X G(r3,r4;w2)G(r4,r2;a;-a;i -1-^2) 

X ^(ri - r4)w(r3 - r2). (26) 

Using a non-interacting Green function of a given spin 



k 



^kn 



(23) 



For a given site R the orbital Green function gna is con- 
fined to the site and it is equivalent to a core state Green 
function. 



GO(r,r':c.)^y^^ii^M(£l + 



n m 

(27) 

we can perform the frequency integral over uj2 using 
Cauchy's theorem by closing the contour either in the 
upper or lower plane: 



C. Theoretical justification of self-screening 
correction 

Here we show that removing the self-screening terms 
in the self-energy is partially equivalent to adding vertex 
corrections in the form of exchange diagrams. It can 
be shown that the self-screening terms are cancelled by 



^G''(r3,r4;w2)G°(r4,r2;w - uji + ^2) 

ZTT 



occ unocc 



^ J>n(r3)<^*(r4)^m(r4)¥'m(l-2) 



w — wi + e„ — £„i -f- iS 
<^rn(r3X,(r4)(^„(r4X(r2) ■ 



UJ — UJi + £r 



iS 



(28) 
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Similarly, integrating over uJi we find 
S"(ri,r2;c.) 

occ unocc „ 

= - ^ ^ / (fir3(fr4^ w(ri - r4)u(r3 - ra) 

^fc (r 1 ) (r3 ) (/?^ (rg ) ip*^ (r4 ) v?„ (r4 ) y?; (ra ) 

c/jfc(ri)v3*(r3)(^„(r3X(r4)(p„(r4X,(r2) 




n m 
occ 



a; - 6fc H- e. 



i(5 



(29) 



The second-order direct self-energy is 



X v{xi — X4)v{x3 — X2), (30) 

and we have considered the direct term with all G having 
the same spin since this is the term that contains self- 
screening. Its Fourier transform is given by 



EVi,r2;-) =|dV3dV4|^|^ X G(r4,r2;a.4) 

X G(r3, r4; uj2)G{r4, r3; w - wi + ^2) 

X ^(ri - r4)u(r3 - r2), (31) 

which can be calculated analytically as in the exchange 
case yielding 



S''(ri,r2;c^) 

OCC unocc » 

^ ^ / Sr^dPr^ v{ri - r4)i;(r3 - V2) 

(y9fc (r 1 ) (y9* (r 2 ) (p„ (r3 ) <^;, (r4 ) V3„ (r4 ) (ra ) 



n ra 
occ 




uj — ek+ Em — e-n — iS 
V'kiri)(pl{r2)(pn{r3)Vni'^4)'Pni{r4)ip*„,{r3) 



UJ - Ek + en - 



iS 



(32) 



Comparison between S'' and reveals that the self- 
screening terms n = in E'' for an occupied cpk, 



fkiri)(pl{r2)(pmir3)ip;^{ri)(pkir4)plir3) 

UJ - ek+ Em - £k- iS 



(33) 



where ipm is unoccupied, are cancelled by the correspond- 
ing terms in S^. Similarly for the case when ipk is unoc- 
cupied. Thus we see that by removing the self-screening 
terms from the direct self-energy we effectively include 
the exchange self-energy. 



III. THE RANDOM-PHASE APPROXIMATION 
WITH SELF-POLARIZATION CORRECTION 

In the previous section we have developed a scheme 
for removing the self-screening in the random-phase ap- 
proximation (RPA)ii in relation to the GW approxima- 
tion. When considering the propagation of an electron or 
a hole that is screened by the surronding electrons, the 
electron or hole in question should not participate in the 
screening process. Here, we apply an analogous idea to 
the case where the perturbation is not due to an electron 
or a hole but to a dipole or an electron-hole excitation. 

In the RPA the polarization is given by 



P(r,r» = ^ 

^Pa(r,r';cj) 



d„(r)d;(r') d„(r')d*(r) 



A, 



(34) 



= Ejn - En ~ i5, £,„ > ^i, En < fJ- (35) 

The index a includes the spin. The response function is 
given by 



R= [l-Pv]-^P 

= P + PvP + PvPvP - 



(36) 



We can think of [1 — Pv\~^ = as a screening factor 
that screens the bare polarization P which consists of 
electron- hole excitations {pa}- We observe that a given 
electron-hole excitation pa generates via the Coulomb 
interaction screening polarizations that include itself be- 
cause P contains Pa- To eliminate this self-polarization 
we therefore calculate the self-polarization corrected re- 
sponse function as follows: 



where 



i? = ^[l-P„i;]- 



Pa = p- Pe- 



rn 



(38) 



Physically this means that a particular polarization p^ 
should not participate again in the screening process so 
that it should be subtracted out from P. To distinguish 
it from self-screening, we have referred to this type of 
process as "self-polarization" although in essence it is 
also a self-screening process. 

Analogous to the self-screening correction described 
before, the self-polarization correction may be regarded 
as an appproximate way of including the exchange di- 
agrams. Consider the first order direct and exchange 
terms. The direct term is given by 
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where 



Pd{xi,X2) = - J d'^X3d'^X4G{x3,xi)G{xi,X3) 

X v{3-A)Gixi,X2)G{x2,Xi). (39) 
For the exchange term we have 

Px{xi,X2) = J d'^X3d'^X4G{x4,Xi)G{x2,X4) 

xGix3,X2)Gixi,x3)vi3-'i). (40) 

Writing the Green functions in Fourier representation 
yields 

Prf(ri,r2;w) 

= - J d^r3d^r4 J ^G(r3, ri; tJi)G(ri, rs; wi + uj) 
-^G(r2, r4; W3)G(r4, r2; W3 + a;)u(r3 - r4) (41) 



X 

and 



v^JM^ J d\d\'^^{r)^*{r)v{r-r'M{r')^i{r'). (45) 

The two self-polarization terms, corresponding to n = n' 
and m = m' in the second and third terms of (j43l) . are 
cancelled by the corresponding terms in P^;. 



IV. APPLICATION TO A MODEL HYDROGEN 
DIMER 

A. The HOMO-LUMO gap in the conventional 
GWA 

Consider a model hydrogen molecule with one orbital 
centered on each atom. The two orbitals centered on 
different hydrogen atoms, ipi and ip2, are normalized but 
not in general orthogonal: {ipi\ip2) ^ 0. The one-particle 
eigenfunctions arc the bonding and anti-bonding states: 



Pi;(ri,r2;a;) 

= J d^r3d^r4 J ^G(r4,ri;wi)G(ri,r3;wi -f-w) 



2tt 



G(r2,r4;tj3)G(r3,r2;a;3 -I- w)t;(r3 - r4). (42) 



Using a non-interacting Green function of a given spin 
yields, using the convention that repeated indices are 
summed and n, n' refer to the occupied orbitals whereas 
TO, to' to the unoccupied orbitals, 



Pd(ri,r2; w) 



ip*Jri)(pmiri)Vnm,n'Tn"Pn'{r2)'P*m' (''2) 
(w - £,„ -I- £„ -I- iS){uJ - Srn' + £«' + iS) 

_ Vn (ri )(/?m(ri )Vn7n,m'n' ^m' i'r2)'P'^' (l'2) 
{uj - Em + e„ + iS){uJ + Sm' - £«' " iS) 

_ (ri ) {ri)Vmn,n'm' fn' (r2 ) 'fil,, (r2 ) 

(CJ + Srn - Sn - iS){uJ - £m' + £n' + i6) 
, '^m(ri)<y5n(ri)Wm„ .m'n' 

{uj + Em - £n - + - £„' - iS) ' 

(43) 

For the exchange term we obtain for a given spin 



Pa;(ri,r2;w) 



V'n' (1-2)9'™' (1-2) 

{lu - e„i + £n + 1(5) (cj - £„i' + £„' + i5) 

(Pm'{'r2)'Pn'(Y2) 
{uj ~ Sm + £n + iS){uJ + Em' - £«' " iS) 
(ri (ri )Vmn'.nm' fn' (r2)t/?^' (r2) 
(CJ + Em - £n - iS){uJ - Em' + £„' -|- iS) 
fmi'^l) fn {ri)Vmm' ,nn' ^m' (r2 ) ^n' {^2 ) 



{uJ + £. 



i6){uj + £„ 



Z(5)' 

(44) 



<P2j 



(46) 



(47) 



with eigenenergies respectively Eb and ea- The indices 
A and _B include the spin function a and /3. These two 
eigenfunctions are orthonormal. We may assume that 
4>i and 4>2 are real. The two electrons occupy the bond- 
ing state with up and down spin. The non-interacting 
Green function (the up and down spin Green functions 
are identical) is given by 

GO(r,r';^) ^ MllMll + MeIME^, (48) 
UJ — E B — io UJ — Ea + io 

where the one-particle Hamiltonian is taken to be the 
Hartree one. The HOMO-LUMO gap in the Hartree ap- 
proximation is 



A" = £^ - £B = 2t, 

where the hopping integral is given by 



t 



1. 



Vil 



Vh\iP2 



(49) 



(50) 



The onsite and intersite Coulomb interactions are respec- 
tively 



Uo ^ {^l\v\^l) = {^l\v\ipl) 
Ui = {^l\v\^l). 



(51) 
(52) 
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{ipiip2\v\ipitf2) and {ip\\v\ipiip2) are neglected since they 
are much smaller compared with Uq and Ui. 

First, let us calculate the exchange contribution: 



^^(v,v') = -v{v-v')^b(j)^b{v'). 



(53) 



The matrix elements in the bonding and anti-bonding 
states are 



(V'bIS'^IVb) = -2(t^o + c/i), 



:(C/o-C/i). 



The HOMO-LUMO gap in the HFA is therefore 



(54) 
(55) 



Using 



Ag 



< 1, 
1 



/(I 



2{Uo-Ui} 
As 

Uo-Ui 
2t ■ 



(64) 



VF'=(r',r;w) = / d^rid^r2v{r' -ri)R{ri,r2;uj)v{r2 - r), 

(65) 

the correlation part of the self-energy E'^ can be calcu- 
lated analytically to yield 



,HF 



2t + Ui. 



(56) 



We now proceed to calculate the correlation part of the 
self-energy. The polarization function can be written in 
the form 



P°{r,r';Lo) = ^B(r)^A(r)P°(c^)^B(r')VA(r'), (57) 
where 



P"(c.) = 2( 1 



1 



with 



e + iS UJ + Ae — iS 



Ae ^ ea- £b- 



(58) 



(59) 



The factor of 2 in ([55)1 is due to the sum over spin. Using 
the RPA equation in and solving it by iteration, it 
is straightforward to see that each term in the iterative 
solution can be written in the same form as P*^ so that 
the response function can also be written as 

R{r,r';aj) = ^PB{r)Mr)R{iu)i,B{r')^Air'). (60) 
R{uj) can be calculated algebraically and it is given by 



Riuj) 



2r 



2r 



uj-AE + i6 oj + AE-i6 



where 



AE= J{Ae)^+4vBA,BAAe, 



(61) 



(62) 



(63) 



J:^aw{r,r';u:)^t J ^G°(r, r'; + a;')W-=(r', r; w') 
Ai(r,r') , A2(r,r') 



oj + AE-EB-iS oj - AE - EA + iS 

(66) 



where 



Ai(r,r') = 2rV'i3(r)?/'B(r') / d^rid^r2v{r - ri) 



X tpB{ri)iJA{ri)'4'A{r2)'4'B{r2)v{r2 -r'), 



A2(r,r') =2rV'A(r)V'A(r') / d^rid\2v{r - r^) 



X ipB{ri)-ipA{ri)ipA{r2)-ipB{r2)v{r2 - r'). 



(67) 



(68) 



(V'B|E^^(cj)|V'i3) 



'^''"^BB.AB 



'^''""\b,ab 



Lj + AE-EB-iS uj-AE-EA+iS 
1 rjUo - Uif 
2uj-AE-EA + iS' 



(69) 



{^A\^^awi^)\4'A} 



'^^^\b.ab 



'^''"'^\a,ab 



uj + AE-EB-iS uj-AE-EA + iS 
1 r({/o-C/i)2 



(70) 



2 cj + AE -EB -iS 

Adding to the Hartree-Fock gap in ^ the HOMO 
LUMO gap in the GWA is therefore 



^GW = 2t + Ui 



rjUo - Uif 
Ae + AE 



(71) 



It is interesting to note that correlation effects increase 
the Hartree-Fock gap, counter to the usual expectation. 
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B. The HOMO-LUMO gap in the GWA with 
self-screening correction 

Let us now apply our GW with self-screening correc- 
tion scheme. For an electron of a given spin in the bond- 
ing or anti-bonding state the screening is provided by 
the other electron with opposite spin, as it should. Thus, 
the polarization is half of the polarization without self- 
screening correction. The calculation proceeds as in the 
previous section and we obtain 



Lj - AE + iS 



AE-iS' 



where 



AE= J{Aey+2vBA^BAAe, 



Ae 

Xe' 



(72) 



(73) 



The correlation part of the self-energy with self-screening 
correction is 



-'GW-SS 



(r,r';a.) = - 



1 Ai(r,r') 



2uj + AE - SB -iS 
I 1 A2(r,rO 
2uj - AE - EA + iS'' 

where Ai and A2 are given in (I57|) and Thus, 



(74) 



2 

_ ''"^BB,AB ^ 

CO + AE — £b — iS 
^ 1 r(t/o - U^f 
Auj - AE - £A + i5' 



(V'aISgh'-ss('^)IV'a) 



2 



- AE - EA + iS 



(75) 



■^AB,AB 



^AA.AB 



AE -EB-iS 
rjUp - Uif 
4a; + AE - eb -iS' 



1 



AE -EA + iS 



(76) 



Taking into account the Hartree-Fock gap in ([55]) . the 
self-screening-corrected GW HOMO-LUMO gap is there- 
fore 



A 



GIV-SS 



2t+Ui + 



rjUo - Uif 
2{Ae + AE)' 



(77) 



where 



C. Exact solution in atomic basis 

We consider configurations with total Sz — 0- In this 
case the Hamiltonian is given by 



-t 



H = 



2eo + Ui 

2eq + Ui t t 

-t t 2eq + Uq 

-t t 2eo + Uo 



which can be solved analytically. Since Eq appears only 
in the diagonal element, we may set it to zero. Choosing 
Eq = 0, the ground-state energy is given by 

EoiN) ^liUo + Ui) - ^^{Uo - Uif + 16^2. (80) 



To calculate the bonding-antibonding or HOMO- 
LUMO gap, we need to consider the ± 1 problems. 
For the one- and three-electron problem there are only 
two configurations. The eigenvalues are 



Ei^2{N +\) = iEo + Un 
Ei^2{N-l)=Eo±t. 



2Ui ± (-i). 
The exact HOMO-LUMO gap with t > is 



(81) 
(82) 



^cxact = Ei{N + 1)- 2Eo{N) + Ei{N - 1) 

= -2t +Ui + v/(C/o - Uif + 16^2. (83) 

It approaches 2t as C/0,1 — > 0, as it should. In the weak 
or moderate coupling regime where {Uq — Ui)/2t < 1 and 
the gap is given by 



^cxact = + C/i + y'iUo - UiY + 16t2 



(84) 



This is the same as the gap in the GW scheme with self- 
screening correction up to order [{Uo — Ui)/2tf: 



A 



GW-SS 



2t + Ui 



2t + Ui 



rjUo - Uif 
2{Ae + AE) 

t f Uo-Ui 
2t 



(85) 



Ae 

'ae 



1 



U^-Ui 
Ae 



-1/2 



(78) 



It is shown below that this is the same as the exact re- 
sult up to second order in (JJq — Ui)/2t in the weak to 
moderate coupling regime where {Uq — Ui)/2t < 1. 



V. CONCLUSION 

We have proposed a scheme for taking into account 
self-screening correction within the GWA. The scheme in- 
troduces orbital and spin dependent screened interaction. 
While this is not necessary in theory, the introduction of 
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orbital and spin dependence W within the GWA captures 
the essential physics better and improves the self-energy 
without resorting to complicated vertex corrections. This 
is analogous to the introduction of the spin variable in 
the spin density functional theory. The scheme is justi- 
fied theoretically by showing that the self-screening tcirms 
are indeed cancelled when exchange diagrams beyond the 
GWA are considered. When applied to a model hydro- 
gen dimcr, the scheme reproduces the exact result in the 
weak to moderate coupling regime. Work is now under 



way to apply the scheme to real systems. 
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